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OO I Abstract 

3^ ' The planar Navicr-Stokcs equation exhibits, in absence of external forces, a trivial asymp- 

*vj [ totics in time. Nevertheless the appearence of coherent structures suggests non-trivial interme- 

diate asymptotics which should be explained in terms of the equation itself. Motivated by the 
separation of the different time scales observed in the dynamics of the Navier-Stokes equation, 
we study the well-posedness and asymptotic behaviour of a constrained equation which neglects 
the variation of the energy and moment of inertia. 
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1 Introduction 



Consider the two-dimensional Euler equation in vorticity form 

{dt + u-V)uj{x,t) = 0, xG]R2 (11) 

where the divergence free velocity field u is given by u = V ■0, ip = — A~"^u;. Explicitely, we can 
write: 

1 X"*" 1 

>i u = K*uj, K{x) =V-^g = - — -r-r^, g{x) = - — log\x\. (1.2) 

m : ' ' 

f^ , The rigorous justification of the formation of coherent structures in two-dimensional fluid-dynamics, 

Cn ■ which is observed in real and numerical experiments (see e.g. [22]) remains a widely open problem. 

An attempt to justify the appearance of these coherent structure is due to Onsager [27], see also 
|17| . |25j . and [TI] for a recent review. The main idea is to replace the incompressible Euler equation 

5o i by the system of A^ point vortices and to study the Statistical Mechanics of these point vortices. 

^D ■ In the mean field limit N — > -|-oo, the Gibbs measure associated to the point vortices concentrates 

to some special stationary solutions of the Euler equation (called mean field solutions), we refer to 
[Z]) [H] ; [E] , [E] for the rigorous justification. These states are under the form: 






yiR2 



In this last expression, Z is a normalization factor to have J lv = 1 and b real and a < are 
parameters. From a mathematical point of view this equation enters in the framework of a general 
class of nonlinear elliptic equations given by 



^bip+a\x 
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AV = -^— (1.4) 
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which has been studied in [7], [8]. Nevertheless, there is no justification of the fact that, among the 
infinite number of stationary stable solutions of the Euler equation, the mean field solutions play 
indeed a special role in the dynamics. 

Another justification of (|1.3p could come from the intermediate asymptotic behaviour of the 
two-dimensional Navier-Stokes equation which, in vorticity form, reads 

{dt + u-V)io{x,t) =uAuj{x,t), x£R^, (1.5) 

where i/ > is the viscosity coefficient. Indeed, due to the dissipation term in the right hand side of 
eq.n (jl.Sp . the asymptotic behaviour of the solutions is trivial, namely, when t — > +oo, uj{x,t) — > 
pointwise and in the L^ sense, for p > 1. Consequently, the states (jl.3p could play a part only in 
the intermediate behaviour of the equation before the dissipation scale. To give a more quantitative 
description of this idea, it is usefull to recall that the solutions of (jl.3p can be studied through 
a variational principle: the radial solutions are obtained as minimizers of the Boltzmann entropy 
(this is proven in the Appendix) 

S{uj) = / ujlogu) dx 
under the constraints 

E{lo) = E, M(u;) = 0, I{uj) = I, u; > 0, ujdx = l, 

for some fixed E and /, where the energy E{uj), the center of vorticity M{uj) and the moment of 
inertia I{u)) are respectively given by 

E{uj) = - il^iodx, M{uj) = / xu>dx, I{u!) = - \x — M\ oj. 

Note that it is always possible to choose the coordinates so that M{lo) = 0. 

It is thus interesting to study how these quantities, which are conserved by the Euler equation 
(jl.ip . evolve under the Navier-Stokes flow. At first, it is well known that the Navier-Stokes equation 
preserves the nonnegativity and that J u> and M are conserved. Consequently, throughout this 
paper, we will focus on non-negative solutions which are normalized such that J to = 1 and M{uj) = 
0. Next, we also observe that 

2 A, ^ f IVWI 



I{uj) = 2v, E = -v \ Lo\ S{u) 



UJ 



It is easy to see that I can be considered as constant for times t « u~^ . Moreover, it is likely that 
in certain cases (see for instance [22]) the energy dissipation rate is much smaller than the entropy 
dissipation rate. Coming back to an attempt to justify eq. (|1.4p in terms of the Navier-Stokes 
evolution, the first naive remark is that, if the energy and the moment of inertia are assumed to 
vary on a long time scale, they can be considered as constant in a first approximation. On such a 
time scale, the motion should be governed by a master equation, which modifies the Navier-Stokes 
equation leaving constant both energy and moment of inertia, but retaining all the other features of 
the Navier-Stokes dynamics. Therefore such a master equation, dissipating the entropy at constant 
energy and moment of inertia, would lead to the solution to eq.n (jl.4p as t — > oo. By using a 
recent geometric gradient flow characterization of the Navier-Stokes equation (see (34j for example) 



connected with the mass transport problem and the associated differential calculus introduced in 
|29j (see also [1]) we have derived such an equation in [10]. In this framework the Navier-Stokes 
equation can be written as a differential equation for a vector field which is the sum of a dissipative 
part, which is the gradient flow of the entropy, and a conservative part, corresponding to the Euler 
equation, which is the orthogonal gradient of the energy. The equation we were looking for was 
obtained by keeping only the orthogonal projection of the vector field in the tangent space of the 
manifold / =const and E =const. Due to the nature of this decomposition, this procedure modifies 
the dissipative part while leaving invariant the conservative part. Thus we have found the equation: 



dfUJ + u ■ Vuj = I'diviVuj — bajVip — aujx) (1-6) 

I |2 

= i>diY[ujV{loguj — bi/j — a——)^, 
where the Lagrange multipliers b and a are given by 



2Ifuj'^ + 2V _ _ ,^_ 2/u;|VVp + l// 



LV 
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and 

V= ujx-'V'il^= / / uj{x,t)u!{y,t) X ■Vg{x — y) dxdy = —-—. (1-8) 
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A way to validate this approach, is to test it in the following simpler and well understood case. 
A special self-similar solution to eq.n ()1.5p is the so called Oseen vortex: 

Note that this is also a solution to the heat equation. It was shown by Gallay and Wayne [12] that 
this solution describes the long time asymptotic of the Navier-Stokes equation in L^. Indeed, with 
the change of variables 

e=^=L=; r = log(l+t), oo{x,t) = {l + t)-'w{^,T), 

\' I + t 

the Navier-Stokes equation in the new variables reads : 

drW + V ■ V^w = uA^w + uV^ ■ ( -S,w ] . (1-9) 

It is possible to show that w ^>- W in L^ as t ^ oo, where W{^) is the rescaled Oseen vortex. As 
a consequence the Oseen vortex can be thought as characterizing an intermediate asymptotics for 
times vt « 1. 

This analysis enters perfectly in the context of the projected gradient flows. Indeed imposing 
the Constance of / in the Navier-Stokes equation we find 

dtUJ + u-Vuj = I'Auj + i^-V ■ (ioxY (1.10) 

that is eq.n (jl.9p for / = 2, as well as the rescaled Oseen vortex VF(^) is a Mean Field solution for 
P = 0. This suggest to impose also the Constance of E in the attempt of outlining what happens 



before the occurrence of the Oseen vortex. Indeed one could argue that / is more robust than E 
in many interesting physical situations. If so eq.n ()1.6p should be more appropriate on the time 
scale when E is practically constant, while eq.n (|1.9p should describe the fluid when E start to be 
dissipated at constant /. 

The aim of this paper is the mathematical study of equations (jl.6p . (jl.7p . For more details on 
the derivation of these equations and of the physical motivations, we refer to [10]. As explained 
in [To], the procedure of constraining a diffusion equation is highly non unique. Nevertheless, an 
interesting feature of eq.n (jl.6p is that it can be obtained by many different methods. As already 
explained, it appears naturally by using the geometric structure of the Navier-Stokes equation. 
Moreover, it was noticed in [10] that eq.n (II. 6p can also be obtained by constraining the stochastic 
vortex dynamics which is a finite dimensional approximation to the Navier-Stokes equation (see 
|18j . |19| [28]...). Indeed (at a formal level), it is shown that the stochastic process for a system 
of N stochastic vortices, once constrained on the E =const manifold, produces, in the mean-field 
limit, exactly eq.n (jl.6p which turns out to be compatible with this particle approximation. 

Equation (II. 6p however is not new. It was previously derived by Chavanis in [5] and [6] following 
a completely different approach based on the kinetic theory of (deterministic) point vortices. 

We finally note that the projection on the manifold / = const according to the gradient notion 
used here, has been considered by Carlen and Gangbo in [3] for a different class of equations. 

A different approach to the one of Onsager [27] in order to understand the coherent structures 
arising in 2D flows was proposed by Robert and Sommeria [31] and Miller [2l]. The equilibrium 
solutions are obtained by using the maximum entropy principle over a state space formed by a 
selected family of possible values of the vorticity. Note that this procedure preserves all the Euler 
invariants so that, as far as the equilibrium is concerned, the Robert-Sommeria-Miller theory is 
quite different from the approaches, as the present one, based on the mean field equation. As 
regards the dynamics, a class of master equations leading to such equilibria has been introduced 
in [32]. One of them, which has some formal similarities with our model, has been systematically 
investigated from a mathematical point of view in [23]. We remark that such an equation exhibits 
a maximum principle leading to useful a priori estimates, for instance the L°° norm of the vorticity 
is uniformly bounded while, in our case, we do not have such a priori control. 

Our aim is to establish global existence results for (II. 6p and to study the asymptotic behaviour 
of global solutions. There are two main difficulties. The first one is that eq.n (jl.6p makes no sense 
whenever the denominator in the definition of a and b (see eq.n (|1.7p ) vanishes. Note that by the 
Cauchy-Schwarz inequality we have 



V^ = ( f uJX-V^y <2I /"w|V^p 



and thus this denominator is always non-negative. Nevertheless, when cvX/ip and ujx are collinear, 
it vanishes. This happens for the one-dimensional family of circular vortex patches: 

1 

where Xb{o,r) is the characteristic function of B{0,R), the disk of center and radius R. Indeed, 
we have: 

The other difficulty is that b is well-defined if lj € L^ but there is no a priori estimates available 
for the L^ norm of the vorticity. The only a priori information we have at our disposal are that E 



and / are conserved (note that in this setting E is not very usefull since the energy has no sign) 
and that the entropy S decays. Indeed, we formaUy have 



dt 



= —u I Loy log uj -yl^log uj — bip — a——) 
f kP |2 



This identity can be checked by direct computation and it is obvious by using the geometric inter- 
pretation of the equation (see [inj ) • Note that this identity is also usefuh to guess that asymptotic 
states should be given by the mean field solutions ()1.3p since the entropy dissipation vanishes 
precisely on these states. 

Because of these difficulties, we will be able to get global existence results only for data suf- 
ficiently close to a mean field solution. Nevertheless, we point out that our smallness constraint 
is independent of the viscosity parameter v. It remains an open problem to establish if eq.n (II. 6p 
can produce a singularity in a finite time and in particular if the L^ norm of lo can blow up. Note 
that if we consider ()1.6p without (jl.7p i.e we consider the equation with some given parameters 
a < and b fixed, it is easy to establish the existence of solution which blow-up. Indeed, since 
the inertial term does not play any part in the estimates, all the result established in [3] for the 
Keller-Segel equation remains true for this equation. In particular we have that the evolution of / 
which is nonnegative is given by 

/ = aI + (2--^) 

and hence, if b is larger than Svr, the solution must blow up in finite time. It would be very 
interesting to know if there is a nonlinear stabilization for (jl.6p , (jl.7p . 

The paper is organized as follows: the global existence proof is presented in Sect. 3 after that 
some preliminary steps are discussed in Sect. 2. Sect. 4 is devoted to the proof of U* estimates 
necessary both for the existence part and the asymptotic behavior discussed in Sect. 5. Finally, 
the Appendix is devoted to the study of equation ()1.3p and the connected variational principles. 
The main ideas are in [7], [8], but the adaptation of these results to the M^ case requires some care. 

2 Preliminaries 

Let us introduce the submanifold of probability densities 

M{E,I) = iuj, uj>0, uj = l,E{uj) = E,I{uj) = l\ (2.1) 

for some fixed E and / > 0. Next, by using Theorems [T2l [T^ in the Appendix, we denote the 
unique minimizer of the entropy functional S{co) on Ai{E, I) by ojmf- Note that ojmf is a (radial) 
solution to eq.n (jl.3p with parameters b and a that we denote by 6mf and umf respectively. 

2.1 A stability property 

We first establish a crucial result asserting the continuity of the L^ norm with respect to variation 
of the entropy, in a neighbourhood of cjmf, in the manifold M{E, I). 



Theorem 1 For any e > there exists 5 > such that, for all uj G M{E, I) for which 

Siio) - SicoMP) < S, (2.2) 

then 

\\uj - ujmfWl^ < £■ (2.3) 

We remark that Thm 1 provides a proof of the Lyapounov stabiHty of lomf with respect to the 
Euler flow (by virtue of the time invariance of S{uj) ). The stability of ujmf was also proved in 
[7] by using the arguments of [20] in which all the conserved quantitities of the Euler equation are 
used. 

Proof of Theorem [Tl Assume, by contradiction, the existence of a sequence uJn G M{E,I) 
and e > such that 

lim5(u;„) = S{u;mf) and ||u;„ — ujmfWl'^ > £• (2-4) 

n 

Thanks to the entropy bound, we can find a probability distribution to G L^ such that, up to 
the extraction of a subsequence, lim„ a;„ = 6<j in the sense of the weak convergence of measures. 
Next, we also have that 

\imE{uJn) =E{uj) =E 

n 

(see the proof of (jA.lip in the Appendix ) and that 

lim/(LL'„) = I > I{uj), limS^iOn) > S{uj) 

n n 

by convexity. Let now umf, ^mf be the multipliers for which ujmf solves eq.n (1.16) with those 
values of parameters, and F(^aMF,bMF) ^^^ ^^^^ nergy functional (see (jA.lh for the definition). We 
get, since omf < 0, that 



^{aMF,bMF)(^) ^ lim \^S{uJn)-bMFE{uJn)-aMFl{uJn) 

= {S{iOMF) - bMFE{cOMF) - CLmfH^Mf)) 
= ^{aMF,bMF)i^MF)- 

Since lomf is the unique minimizer of F{bMF,a.MF) (^^^ Theorem [T2|) . it follows that lo = ujmf- As 
a consequence, we also get that 

\\mI{uJn)=I{uj)=I{uMF)- (2.5) 

n 

Finally, we consider the relative entropy 

S{uJn\u:MF) = /w„log( ^) 

J ^MF 

= S{u)n) - S{uJmf) + bMF / i^MF " (^n) IpMF + ^MF U (^mf) - I{^)) ■ 

Now, we observe that S{uJn) — S{lu) goes to zero thanks to (]2.4p and that 

bMF / y^MF - ^njIpMF + dMF \^I{^Mf) " H^)) 



also goes to zero by weak convergence and (j2.5p . Consequently, the relative entropy 5(u;„|u;m_f) 
goes to zero. Thus we conclude by using the Csiszar-Kullback inequality 

||w„, — iOMpWii < 25(u;„|a;MF) -^ 
which yields the desired contradiction. ■ 

2.2 Properties of the coefficients a{uj),b{uj). 

Our next step is the study of the properties of a and b. We shall use the notation 



2//w2 + 2y 2/w|VVp + l^/ 



UJ 
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where 



L»(w) = 2/ /"a;|VV'p-V^^ V = - — . (2.7) 
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As we have seen in the introduction, one of the main difficulty is that the denominator D{lj) may 
vanish for a vortex patch. 

Before stating the result, we shall recall a usefull set of inequalities which will be used throughout 
the paper (see e.g. [33]...). 

Lemma 2 (Useful inequalities in M^) 

i) Sobolev-Gagliardo-Niremberg : The following inequalities hold for some C > 0: 

I |u;| 1^2 < Clival lii (2.8) 

\\uj\\l2 <C\\uj\\li\\Vuj\\l2. (2.9) 

ii) Biot et Savart law : Let u = K * uj with K defined by (II. 2p . then we have 

l<p<2, 2 <q < +00, - = -, ||n||L9 < CIIwIIlp, (2.10) 

q p 2 

l<p<2, 2<g<oo, ^ = - + ^—^, \\u\\l^ < CIMIp ||a;||i;", (2.11) 

2 p q 

1 <p< +00, IIVuIIlp < C||w||lp. (2.12) 

Hi) Interpolation in L^ spaces: 

1 — - 
l<p<r <q<+oo, llwllir < 2||(j||2p||w||^-", a = --—^ (2.13) 

We shall prove the following result: 

Theorem 3 Suppose lo £ L^ n M{E, I) for some p > 2. 

1. We have: 

bitOMp) = bMF, a{ujMF) = dMF (2.14) 

and D{ujmf) > 0. 



2. D{uj), a{uj) and 6(u;) are continuous on L^ 

3. Assume that p > 2, then if S{uj) — S{ujmf) is sufficiently small, we have 

\D{oo) - D{ujmf)\ + \a{uj) - a{ujMF)\ + \b{uj) - Hujmf)] < C\\uj - ujmfWIi, (2.15) 

for some a < 1. Here C depends on E,I and \\uj\\lp. 

Proof of Theorem [3l 

We first prove 1. Since ujmf is a solution of the mean field equation (II. 3p . we have 



VuJmF = buF^MF^i^MF + ImF^MFX. 

Taking the scalar product of this equation by 'VipMF and x we find 

bMF / WAfFlVV'A/Fp + clmfV = / wfj^, bMFV + 2aMFl = -2. 
The resolution of this two by two linear system precisely gives that 

b{uJMF) = bMF, a{uJMF) = a-MF- 

Consequently, since the numerator in the definition of b{uJMF) and a{uJMF) is finite, we find that 

D{umf) > 0. (2.16) 

Next, we shall estimate the differences 

ju^-julj,, (2.17) 

/c^lVVp - j UJMF\y^MF?- (2.18) 

By Cauchy-Schwarz, we obtain 

1 (12.17^ 1 < \\uj -ujmfWl'^ (II^IIl2 + W^mfWl'^) 

and hence we find 

IdlUDI <C||w-a;MF|l2i (2.19) 

for some q > 0, by using (I2.13|) with r = 2 and q = I. Next, we split (J2.18P as 

IdHlDI < f{uJ-u;MF)\V^\^+ fiOMFy{lp-^MF)-yilp + 1pMF) 
^ II II II Il2 

^W^MfWl'^ (Iklli^ + ||tiAfF||L4)||ti -UmfWl'^ ■ 

4 

We notice that thanks to (j2.10p . the L norm of the velocity is bounded in term of the Ls norm of 
the vorticity. Therefore, by a new use of ()2.13p . we find 

IdmDI < C (\\lO - COmfWl^ + \\UJ - COmfW^^) <C\\uj-UJMF\\tl- (2.20) 



Next, by using (l2J9l) . (l2:20l) . we get that 

\Diuj) - DiuJMF)\ <C\\uj-ujmf\\'Ii. 
Consequently, we can use Theorem [T] to get that 

D{io) > D{ujmf) - \D{uj) - D{umf)\ > ^I?(wmf) (2.21) 

provided S{uj) - S{lomf) is sufficiently small. Finally, by using (|2:2TD . (12391) . (12:20]) and (I236D . we 
easily conclude the proof. ■ 



3 Global Existence and uniqueness 

We start with a brief explanation about the construction of a classical local solution. Let wq G L^ 
with p > 2 be the initial condition such that ujq G M{E,I). Let comf the unique Mean-Field 
solution associated to Ai{E, I) as above. We assume that S{ujq) — S{u]\fp) is small enough so that, 
by virtue of Theorems [H El we have 

\D{loo) - D{ujmf)\ < -D{ljmf) 

and also 

|a(cJo) - a{u:MF)\ < -\a{uJMF)\, \H^q) - b{uJMF)\ < -|6(a;MF)|. 

This implies that we have the upper bound 

|a(u;o)| + \b{uJo)\ < 2{\a{uJMF)\ + \Kujmf)\) 

and also that 

D{uJo) > -D{ujmf) > 0. 

Note that the positivity of D is important in order to stay away from the singularity. 

For every p > 2, by using a standard iterative scheme, we can easily establish, a local existence 
and uniqueness result for a classical solution uj G C([0, T],LP) n L°°([0, T], L^((l + \x\'^)dx), for 
which 

II'^WIIlp < 2(||a;o||Lp + \\u:mf\\lp), (3.1) 

and 

D{uj{t)) > ^D(wmf), Vt G [0,T]. (3.2) 

Moreover, we can continue the solution as long as the U" norm of lo remains finite and the denom- 
inator D{uj) remains positive. 

Note that L?' seems the natural space for our equation in order to have a and h well-defined. 
Note also that the condition (j3.2p allows to avoid the singularity of D[lo) by Theorem [3l 

Let T > be the maximal time for which, the estimates (j3.ip . (j3.2|) are verified. Our purpose 
is to prove, by a priori estimates, that T = -|-oo. 

To do this we shall use the L^ estimates given by the following Theorem which will be proven 
in the next section 

9 



Theorem 4 Consider a local solution as above of (jl.6p such that 

\a{uj{t))\ + \b{uj{t))\ < Co, Vt G [0,r]. (3.3) 

Assume also that ujq G L^, p G [2, +oo) is a probability density. Then there exists Cp which depends 
only on ujq, Cq and p (and hence does not depend on v and T if Cq does not) such that 

Mt)\\LP<Cp, vtG[o,r]. (3.4) 

We are now in position to get a global existence result by showing that T = +oo. Indeed by 
the H-Theorem (see (1.22)), we have 

S{u;{t)) - S{lomf) < Siioo) - S{u;mf), Vi G [0, T] (3.5) 

and thus by Theorem [H we get that ||a;(t)) — lomfWl^ < e, Vt G [0, T] ( with e independent of T), 
provided that 5(a;o) — S{ujmf) is sufficiently small. Next, we can use Theorem [H Indeed, because 
of (13. 1|) . (13. 2p we have the bound |a| + |6| < Co on [0, T] for some Cq > This yields a control of 
the LP norm of uj with p > 2 , depending only on Cq. In particular, thanks to (j3.2p . we get that 

\D{Lo{t)) - D{u;mf)\ < Ce", \\Lo{t) - lomfWlp < Ce", Vt G [0,T] 
where C depends on Cq only. Consequently, we can choose e sufficiently small to have 

||a;(t)||LP < 2(||a;o||LP + \\ujmf\\lp), D{Lo{t)) > -D{ujmf), Vt G [0,T]. 

and hence T = +oo. ■ 

Therefore we have proven the following global existence result: 

Theorem 5 There exists 6q > (independent of u > 0) such that, for any initial datum loq G 
LP n M{E, I) with p > 2, close to ojmf in the sense that 

S{uJo) - S{ujmf) < 5o, (3.6) 

there exists a unique classical solution 

ojit)£Ci[0,+oo[,LP)nL°°{[0,T],M{E,I)) 

to eq.n (jl.6p with initial datum loq. 

Moreover, we have the Lyapounov stability of ujmf, namely, for any e > 0, there exists S, 
< 6 < 6o such that if S{uJo) — S{ujmf) < S, then 

\\u!{t) — ojmfWl'^ < £, Vt > 0. 

As noticed after Theorem [H lomf is also Lyapounov stable as a stationary solution of the Euler 
equation 11.11 Consequently, we also have the following global stability result between the flows of 
(jl.6|) and (jl.ip in the vicinity of ujmf'- 

Corollary 6 For every e > 0, there exists (5 > such that if uq G L°° nAi{E,I) verifies 

S{uJo) - S{ujmf) < S 

then the global solution uj of Euler equation \l.l\ and the global solution lo'^ of ()1.6p for u > with 
the same initial datum ujq satisfy 

||cj^(t)-w-^(t)||ii <e, Vt>0. 

Note that this global approximation property of the Euler evolution (even for z/ large) is of course 
false for the Navier-Stokes evolution. 

10 



4 Propagation of L^ regularity 

In this section we prove Theorem HI 

Our stategy in proving Theorem U] wih be based on weighted energy estimates because in this 
way we can use the fact that the inertial term u ■ Vlo does not contribute. This is crucial in order 
to find estimates independent of u. 

We first focus on the L^ estimate. We shah use the notation 

B = supt^[o,T]iMt)\ + \bm. 
The standard L^ energy estimate for (jl.6p gives: 

|(^IKi)lli2)+H|V^(i)lli2 = z.(6(t)|..(tf-a(t)|u;(t)2). (4.1) 

Next, as in [3], [H], we can use the Sobolev inequahty (|2.9p to get (we recah that J u = 1) 

fio^ = f{Jf <^c(^f\Vio\io-^y < ^CWVioWl^ 

where C is the best constant in the Gaghardo-Nirenberg-Sobolev inequahty. Consequently we get 

|(^ll'^WIli2)+H|Va;(t)||i, <^CS||Vu;(t)||i. + ^5||a;(i)||i, 

where C is an explicit harmless number. Now, let us assume for the moment that CB is sufficiently 
small (less than 1/2 for example), then we can deduce that 

^(^IkWIli^) + ^||V.;(t)||i. < uBMrnl,. (4.2) 

If we directly integrate this differential inequality, we still cannot conclude. Indeed we shall find 
that ||a;(t)||2,2 grows exponentially in time and this does not allow us to get a uniform in time 
estimate. Note however that the bad term in the right hand-side of (4.7) comes from the linear 
term V ■ (xu)) in eq.n (|1.6|) . The explanation for the bad behaviour we get is simple: the semigroup 
generated by the linear Fokker-Planck operator 

Llu = AiO + ^V ■ (XLO), 

with i? > 0, is not uniformly bounded in time as an operator in C{Lp'). Nevertheless, it is bounded 
as an operator in C[L^ HL'^, L^). A very simple way to see this property is to use a weighted energy 
estimate. Indeed, multiplying (j4.ip by e'^^*, we find 

^(^ll^Wlli^) + ^-||Va;(t)||i. < -uBe-'^'Mml, Vt G [0,T]. (4.3) 

Now, we can use (j2.9p and the Young inequality to get, for some harmless explicit number C 
(independent of z/) which changes from line to line, 

^iyBe'''^'\\uj{t)\\l2 < C uB e'''^'\\Vuj{t)\\L2 < ^ue'''^'\\Vu{t)\\l2 + C uB'^ e"'^'. 
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Consequently, we can plug this last inequality in (j4.3p to get 

The integration finally gives 

Mt)\\l, <||u;o||i2+CS, VtG[0,r]. 

We now remove the assumption on the smallness of B and prove also the propagation of the L^ 
regularity in the general case p > 2. 

The starting point is to use the idea of [H], [3] in the study of the Keller-Segel equation. For 
K > 1, a parameter which will be fixed later, we define 



niKit) = / {io{t) - K)+dx. 
We note that 

rriKit) < f ioit) < -\ f uj{t)\og{uj{t))dx < -^ / w(t)| log w(t)|. (4.4) 

JLo{t)>K ^ogK J logK J 

Now we can use the following useful inequality : 

Lemma 7 There exists C > 0, such that, for all probability distribution uj, we have 

f u;\ log uj\<S{u;) + C {I + I (u;)) (4.5) 

We postpone the proof of the lemma to the end of the section. 

Thanks to (j4.5p . we find that J uj\ log6<j| is bounded in terms of the initial datum because the 
entropy is decreasing and I{uj) is constant. Thus 

rriKit) < -^Citoo), Vt > 0. (4.6) 

log A 



Next, we can perform a modified L^ energy estimate for the solution of (jl.6p . After a few 
integrations by parts, we find 

I (i Jiu - K)l) + i^ip - 1) Jiu - K)l~'\Viu; - KU\' (4.7) 

= uB(^ j{oj - K)l+^ + {2K +-) + (K^ + 2K) f (u; - Kf-^'^ 



To estimate the first term in the right hand side of (j4.7p . we use the Sobolev-Gagliardo-Niremberg 
inequality (12. Sp and Cauchy-Schwarz. We have 



< C7^^±^m^(t)|(.;-AX'|V(o^-K)H-P 
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This yields, thanks to our assumption (j3.3p and (j4.6p . 

AlogK 
Hence, by choosing K such that 

CoSo{p + lf ^ 1 
41ogK 2 

we obtain 



. B / (c - K)l^' < V ^°^°/^+^^' /(c - iTf -^|V(a; - ir)+|2. 



(p-1) 



vB jiuj- Kf^' < '^^^ 1(0; - K)l-'\V{u - K)^]'. (4.8) 

Note that K depends only on Co, ^o, p and is diverging with p. 
To estimate the last term in (I4.7p . we write 

J JK<u;<K+1 Ju>K+l 

< 1+ [{to- K)l (4.9) 

By plugging (USD and (USD in ([12]), we find 



< vCoiK"^ + 4:K + - /(cj - i<:)!^ + i^2 ^ 2K 
<CoCun\LO-K)l + l 



where, from now on, C is a harmless number which depends only on K and p. Again, we note 
that we cannot directly conclude by using the Gronwall Lemma in the last differential inequality 
because it gives an estimate which is not uniform in time. We now use the technique that we have 
explained in the beginning. We find 



- [e^'- j{u:- K)l) + ^^^e^' y (^ - A')riV(a. - K)+|^ (4.10) 

<C{l + Co)ue'''(^f{oo-K)l + iy 
Next, we can use the inequality (12. 9p to get 

l{co-K)l = l[{u;-K)iy<cl{u-K)l[J\u;-K)l-'\V{u;-K)+\'f. (4.11) 

By using the interpolation inequality (j2.13p of L^" between L^ and L^, we have 

/(w -K)1<C(^ f{io - Kf ^ ^^ 
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and hence, we deduce from (j4.1ip that 

J\co-K)l<c[J\u:-K)l-'\V{co-KU\^)' 

where q is such that p~^ + q~^ = 1. Thanks to the inequahty 

ab<p-^aP + q~^a'i, a > 0, 6 > 0, (4.12) 

we finally obtain 

C(l + Co) jiu; -K)l<^ jioj - K)l-'\V{u; - K^l' + Cp, 

where Cp will now stand for a number which depends only on ujq, Cq and p. By using this last 
inequality in ()4.10p . we finally arrive to 



The integration gives 

[{u; - K)l < Cp. (4.13) 

Now we can conclude as in |14] . [3]. By using the inequality 



x^ < (j^) {x - ly 



for every x > A > 1, we find 



[ujP < f u;P+ f 

J Juj<K Jw 



UjP 



>K 



JK<uj<\K Juj>XK 

< KP-' + (XKy-' + ^^y~' J (u - K)l. 

This ends the proof of Theorem [J] ■. 

It remains to prove Lemma [7] which is a classical estimate we present for completeness. 
Define uJ = a;l||(^|<i} Since we have 

/ u;|logLL)| = ^(lj) - 2 / ZJloguJ, (4.14) 

it suffices to find a bound from below of JuJloguJ. By using the fact that the relative entropy 
between two probability measures is non negative (this is an easy consequence of the Jenssen 
inequality) we get 

/ (ZJ/m) log ( — 
Jr2 ^J_, 



27r 
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e 2 



> 



where m = J^2 uj < 1. Then we get 

/ UJ log io > —mliuj) + m log m log m > —lioj) + log 

J 27r 27r e 

and hence we get (|4.5|) by using this last estimate and (|4.14p ■ 

5 Asymptotic behaviour 

In this section we investigate the asymptotic behaviour of the global solutions given by Theorem 
m More generally, one can consider a global solution uj{t) of (jl.6p such that u> G C([0,+oo[, L^ n 
L-^((l + |a;p)dx) and such that uj{t) G A4{E,I) and which satisfies for some C > 0, the uniform 
estimates 

\\io{t)\\L2 < C, \b{uj{t))\ + \a{uj{t))\<C, Vt>0. (5.1) 

The main result of this section is given by the following Theorem. 

Theorem 8 Let uj{t) a global solution of (II. 6p as above which satisfies (15. ip . Then uj{t) converges 
in L^ , ast ^ oo, to the unique solution ujmf £ M{E, I) of the associated microcanonical variational 
problem. 

Note that the solutions constructed in Theorem [5] satisfy the estimate (|5.ip and hence their 
asymptotic behaviour is given by Theorem [HI 

Proof of Theorem [8] 

The first step consists in proving that the orbit {Lo{t)}t>o is relatively compact in L^ and uniformly 
bounded in L°°. Before, we need to study the evolution operator generated by the non-autonomous 
Fokker-Planck type operator 

L^u> = V ( Aw + 7(t) V • {xijS) 



where 7(t) is a given continuous curve. Denoting by 5'^(i,r)tJo the solution of 

dtuj = L^w, t > T, w(r) = Wo, 
we have the following estimates 

Lemma 9 Suppose that for all t > 0, \^{t)\ < Kq, for some Kq > 0. Then, there exists C > 
independent of u > and such that for [p, q, r) G [1, +00]'^ we have: 

\\S^{t,T)uj\\LP <CKq^~ — '- -^\\uj\\li, - + - = 1 + -, (5.2) 

\\VS^{t,T)uj\\LP<CK^ ^ — -^^||a;||L,, - + - = ! + - (5.3) 

3_i g^ii-o(l-|)(i-r) ^ ^ ^ 

\\S^{t,T)Vuj\\LP <CK^ ^ ^^II^IIl,, - + - = ! + - (5.4) 

for all p £ [1, +00] 
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Proof of Lemma [9] 

A simple computation in Fourier space allows to find the explicit representation 

\x-y? 

5,(t,r)u;(x) =e2'^(^W-^M) ' J ^ ,^, ^ _, ^ u:o{e'^^^^'^-^^^^\) dy (5.5) 



A^v Jl eMB{s)-B{t)) ds 
where B{t) = L 'y{s)ds. The result of Lemma [9] then follows by standard convolution estimates. 



We come back to the proof of Theorem[8l We shall prove that 1 1 Vti;(i)| |j;^i is uniformly bounded. 
We use the same idea as in [12]. Note that the solution of ()1.6p can be written as 



(jj{t) = Sa{t,0)uJo + / SaitjT)^ ■ ( — uuj + ubu w j (t) dr. 



(5.6) 



Moreover, thanks to (j5.ip . we have a uniform estimate on ||a;(t)||^2 and on \a{t)\ + \b{t)\ for all 
times : there exists Cq > such that 

\\Lo{t)\\L2<Co, Ht)\ + \b{t)\<Ko,yt>0. (5.7) 

Consequently, in ()5.6p . we can consider a and b as known and we can use the estimates of Lemma 
m Let us define F{u)) as the right-hand side of (|5.6p . We have 

WFicomiLoo < C(^ - + Co{l + iy) jWucoWL^ds), 

^ a^{t)'i Jo a^{t - s)i ^ 

where a^if) = 1 — e~^'^^o*. Since by (j2.10p . (j2.13p and the uniform L^ bound, we have : 

||^*<-^||l4 < II'^IIl"" II^IIl^ ^ C'II'^IIl°° II'^II 4 ^ C'Collwllioo, 

Li 3 

we finally get 

fffpVKQt rt uKo{t-s) , 

mu;{t))\\L^<C(^—^ + Co{l + u) / — -^Ms)\\L^ds). 

^ au{t)2 Jo au{t - s)i ' 

Consequently, we can set 



[o,t; 
to get 



z(T) = supfe-^^°V(i)^ll^(i)l|j 



z(T)<CCo(l + i^)(l + a,(r)l( / — ,. ds\z{T) 

Next, we notice that 



limauiT)^ f ^3 r=0, 

^-0 Jo a^{T - s)-^ a^{s)2 
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therefore, there exists T(z^, Co) > such that 

z{T{v, Co)) < CCo(l + u) + \z{T{v, Co)) 

and hence, we get that 

||^(i)||^^<£(^^, VtG[0,T(i.,Co)]. (5.8) 

ay{t)-2 



Next, since to estabhsh (I5.8P we have only used (15. 7p . we can consider for every n G N, the solution 
uj of (jl.6p with initial value oj{nT{u,CQ)/2). By the above argument, we get that to satisfies the 
estimate (|5.9p . By uniqueness, we have 



Cj{t) = uj{t + nT{i^, Co)/2), Vt G [0, T(z>, Co)] 

and hence, 

II^WIIl- <C(i/,Co)fl + ^^y Vt>0. (5.9) 

^ a^(t)2^ 

for some C(i/, Co). 

In a similar way we have by Duhamel formula 

Vuj{t) = VS{t,0)+uJo VS{t-T)(-u-Vuj + i^bu^-Vuj + ubu;^]{^)ds (5.10) 

\\VF{uj)\\li<c(-^ + {1 + u)Co f ^ ^ A \\uVu:{s)\\L^+C,)ds 
^ay{t)2 Jo a^t - s)2 

and since we have by (|2.1ip 

\\uVuj\\]^i < C||u||i,oo||Vti;||^i < ClltLiJI^oo ||Vw||^i, 
we get, thanks to (|5.9p . 

||VF(^)||^i < C(-^ + Co / --^—-T ds + C{u, Co) / — ^^— -T l|Vu;(5)Li ds 

^aiy[t)2 Jo a„[t — s)2 Jo a,y{t — s)2ay{T)* 

Consequently, by using the same method as before, we can easily obtain 



\\VLo{t)\\Li <C(i/,Co)(l + ^^),Vt>0 (5.11) 



for some C(z/, Co). 



We now consider 0, the omega limit set of the trajectory {uj(t))t>o- We deduce from the previous 
estimates that the positive orbit {ui{t)}t>o is relatively compact m X = L^{{1 + |2;|°)(i2;) n L^ for 
a < 2. Indeed, since uj{t) G C([0, +oo[, X), it suffices to prove that {uj{t)}t>i is relatively compact. 
The compactness in L^{1 + |a^|") follows immediately from the Riesz-Frechet-Kolmogorov criterion: 
Lo{t) is uniformly bounded in L^, the uniform (for t > 1) bound (j5.1ip gives the equi-integrability 
and we have a uniform bound on the moment of inertia to control the mass far away. Next, thanks 
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to the uniform L°° estimate for t > 1 given by (j5.9p and the relative compactness in L^, we also 
get that {u!{t)}t>o is relatively compact in U^ for every p < +00. 

By the relative compactness properties that we have just proven, we get that 0, is non empty 
and actually made by smooth L^ functions thanks to the smoothing effect of the parabolicity. 
Moreover, we also have that the elements of Q are probability densities. Also, if u; E 0, since there 
exists an increasing sequence t„ such that w(t„,) tends to u; in X, we also have 

E{uj) = lim^(u;(t„)) = E, I{uj) < lim/(tj(t„)) = /. (5.12) 

71 n 

The first equality is proven in the Appendix, see (IA.11|) . Finally, we notice that the entropy 
S is constant on 0. Indeed, if wi, u;2 G ^, we can construct an increasing sequence t^ such 
that uj{t2n) tends to wi and uj{t2n+i) tends to UJ2 almost everywhere and such that there exists 
51, 92 G iH(l + \x\^)dx) n L^ with 

^{t2n) < gi, W(i2n+l) < 52- 

By using that 

uj\ logwl < C(uj'^ + \uj\i) < Ciio^ + (1 + \x\)^ + 



(l + l^l: 

we find by Lebesgue Theorem that 

S{loi) = Iim5'(u;(t2n)), 'S'(a;2) = limS'(u;(t2n+i))- 

n n 

But, since the entropy is decreasing, we also have S{uj{t2n)) > S{uj{t2n+i)) > S{uj{t2n+2)) so that 
passing to the limit, we get 

5(a;i) > S{U2) > S{uJi) 

and hence S{u}i) = S{ijJ2)- 

Finally, we can prove that the elements of Vt are solutions of the mean field equation. If u; G VL, 
consider uj[t) the solution of (jl.6p with initial value lo. By the strong parabolic principle, we have 
that oj{t) is smooth and strictly positive for t > 0. Since $7 is invariant and S is constant on it, the 
entropy dissipation identity (1.22) gives that for t > 

v(loga;(t)-6V(i)-a^)) =0. 
By continuity in time we get that oo actually solves the mean field equation 

u = ^e^^^^-4- (5.13) 

in R^. Finally, by using the result of |26j and Lemma 4.3 of [3], we get that w is radially symmetric. 
Note that, we also necessarily have that a < and 6 < Svr. 

To summarize, we have proven that the omega limit set fi of {u;(t)}t>o is made by probability 
densities which are radially symmetric solutions of the mean field equation (j5.13p with finite energy 
equals to E and finite moment of inertia. Since the entropy separates the radial mean field solutions 
(see Remark [15] in the Appendix), we conclude that Vl consists in a single point. ■ 
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A The Mean-Field Equation and related variational problems 

In this Appendix we collect some useful facts concerning the Mean-Field Equation (MFE) in M?. 
The main ideas are in [7], [8], here, we adapt the results to the R^ case. For the microcanonical 
problem, the strategy of the proof is slightly different, we do not prove directly the existence of 
a solution. We focus on the negative temperature case (which corresponds to 6 > 0) which is the 
most interesting case. 

Definition 10 (Canonical Variational Principle) For a < 0, b > 0, consider the free-energy 
functional 

Fa,b{io) = S{lu) - hE{uj) - al{uj) (A.l) 

defined on the space F of probability densities on M^ for which E{u}), I{uj) and S{uj) are finite. We 
set 

F{a,b) = mi Fa,b{oo) (A.2) 

Definition 11 (Microcanonical Variational Principle) For £■ G M and I > let us define 

M{EJ) = {L0€r: E{u) = E, I{uj) = I] . 

We set 

S= inf S{u;) (A.3) 

The main results of this Appendix are the two following theorems. For the positive temperature 
case (i.e. 6 < 0), F^^a is a convex functional so that there is a unique minimizer. Moreover, there 
exists a unique solution to the MFE [l^ and all the following results are obvious. 



Theorem 12 (Canonical Variational Principle) For a < 0, and < b < 8tt : 

i) there exists (j € F such that Fa^ti^) = F{a,b). Moreover lo is radially symmetric and solves 
the mean field equation (jl.4p . 

ii) There is only one radially symmetric solution of the mean field equation (II. 4p . 

Hi) As a consequence, there exists a unique minimizer lj a, b of F{a,b) overT. 

Remark 13 It is easy to prove that when b -^ Svr, the solutions to the MFE concentrates at the 
origin. Indeed by multiplying the equation by x ■ Vip and integrating by parts, we arrive to the 
identity (that is the same argument leading to the Pohozhaev inequality): 

1-^/=^. (A.4) 

b b 

Hence when b -^ Svr, / — > and the concentration takes place. For b > Stt we do not have solutions. 

As a consequence, we can solve the microcanonical variational principle. 
Theorem 14 (Microcanonical Variational Principle) For a < 0, and < b < 8tt : 

i) F(a, b) is a concave smooth function. 
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a) 

dF dF 

— = -I{u;a,b) < 0, ^ = -E{iJa,b) (A.5) 

oa GO 

in) For I > 0,E eR let us define S*{I, E) as 

S* (/, E) = sup f F(a, b)+bE + al) , (A.6) 

a,b ^ ' 

Denote by a{I , E) , b{I , E) the unique maximizer of ()A.6p . then for any / > 0, i? € M, 
S{I,E) = S*{I,E), and hence S is a smooth convex function. Moreover, the microcanonical 
variational principle admits a unique minimizer Cji^e in ^i,E- Finally Cji^e = ^a{i,E),b{i,E) 
(equivalence of the ensembles). 

Remark 15 

We finally underline that the function I h^ S{E, I) is strictly decreasing (dS/dl = a < 0) so 
that different radial solutions of the MFE with the same energy cannot have the same entropy. 

Proof of Theorem 1121 

By the logarithmic Hardy-Littlewood-Sobolev inequahty (see [2], [9]), we have 

5(u;)-87rS(u;) > -(1 + log^). (A.7) 

Note that we also have the inequality (see (lA.OP ) 

E{u:)>-j-\og{U{u)). (A.8) 

Indeed, we can write 

E = -— \og\x- y\'^iv{x)uj{y) > -— log / \x - y\'^uj{x)uj{y) (A. 9) 

and hence 

E>-^log(4I-2(^J xioY^ >-i-log(4/). 

Consequently, thanks to ()A.7p . (jA.Sp . we get that 

F(,,6)(u;) > -^^ log(4/) - a/ - (1 + logvr) (A.IO) 

and hence, we find that -^"((,,0) is bounded from below. 

Let uJn be a minimizing sequence in T. Up to the extraction of a subsequence, aj„ converges 
in the sense of weak convergence of measures. Moreover, thanks to (jA.lOp . we get that I{uJn) is 
uniformly bounded and, by using again (|A.7p . we also have 

(i_|_)SK,<F(.„) + i:^ 

and hence S is bounded from above. Therefore the uniform integrability given by the bounds on 
5 and / (which yields a bound on Ju;|logu;| thanks to Lemma [7]) implies that u;„ converges to 
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a nonegative function oo. Moreover, the uniform estimate on the moment of inertia provides the 
tightness of the sequence u;„ so that we obtain 



hm uJn = uj 

i.e. oj £T. Next, by lower semi-continuity, we have 

S{lo) < liuiS{u)n), I{oj) < lim I (uJn) 

n n 

and we claim that 

E{uj)=\n-iiE{u:n). (A.ll) 

n 

This proves that F(^a,b){^) ^ hm„ F(a,fe)(u;„) and hence that uj is a minimizer. It remains to prove 
([AH]) . We write 

E{uJn) = -— / log |x - y\iVn{x)uJn{y) dxdy = I{e) + J(e) 
where 
/(e) = --— / \og\x-y\u:n{x)u:n{y)dxdy, J{e) = --— log \x - y\uJn{x)uJn{y) dxdy 

4vr J\x-y\<e 47r J\x-y\>e 

for every e G (0,1). By splitting the integration domain in {iOn{x)ijJn{y) < \x — y\^^} and its 
complementary, we easily get that 

lie) < -C \x -y\~^ log \x-y\ dxdy + 2C( uJn\'^oguJn\) sup Un{y) dy 

J\x-y\<e ^J ' -i: J\x~y\<e 

J\x-y\<e J\x-y\<e 

For the last line, we have used that the entropy and the moment of inertia are uniformly bounded 
in n and thanks to Lemma[71 we also have that /u;„| loga;„| is uniformly bounded. This yields the 
uniform integrability : 

lim swpx / ^n{v) dy = 0. 

"^^0 J\x-y\<e 

Consequently, we get that lim£__»o /(e) = uniformly in n. Finally, by weak convergence, we have 
that 

limJ(e) = — -— / log\x — y\uj{x)uo{y) dxdy 



471" J\x~y\>e 

and since 

1 f 

lim—-— / log\x — y\oj{x)uj{y) dxdy = E{ijj)^ 

e^O A-K J\x-y\>e 

the conclusion follows easily. 

By symmetrizing uj (around the origin) we find that -F(a,fe) is decreasing. Indeed S and / are 
unchanged and E is increasing. Thus uj must be radially symmetric. It is not difficult to show that 
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a; > (otherwise one could find a better distribution as regards the minimization problem). Hence 
uj satisfies the MFE. 

Next we show that such a solution is also unique among all the radial solution to the MFE. 
Setting r = |x| and ilj{r) = ip{x) (by an obvious notational abuse), we have 

i(rV^')' = -e^'/'+t^'. (A.12) 

r 

We are assuming that Z = 2tt J„°° drre^~^^^ = 1, adding, if necessary, a constant to '(/'■ After the 
change of variable t = logr, setting H = bifj + 2t we readily arrive to the following equation: 

H = -F{t)e^ (A.13) 

where 

F{t) = be''—. (A.14) 

We are looking for smooth solutions to eq.n (18) and hence 

lim H = 2 (A. 15) 

t— ►— oo 

as a consequence of the fact that \inir^Q rip' (r) — > 0. H{t) behaves as 2i + x as i — > — oo and x 
must be choosen in such a way that Z = 1. In the new variables: 

97r f'°° 27r 

Z = — dte^F{t) = —{H{-oo)-H{oo)) (A.16) 

b J-oo b 

It is convenient to change the time variable by setting 2t — > 2i — Xi so that the problem can be 
reformulated as 

H = -F{t-^)e^ (A.17) 

H{-oo) = 2, H{t) PS 2t for t -^ -oo. 

Z{x)^0 for x^-oo (A.18) 



Note that 
and 



Z(x)^y for x^+oo. (A.19) 

Eq.n (jXj) is obvious, while eq.n ()A.19p comes out by integrating the Hamiltonian system 

H = -e^ (A.20) 

for which, the energy conservation yields H{oo) = —2. 

Since ^ > 1 the value Z = 1 is certainly taken, at least once. In order to get uniqueness it 
remain to show that x -^ ^ix) is a monotone function, actually it is not decreasing. 

Defining 

G = H+-e^^'-^\ (A.21) 
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we find the following set of non-autonomous equations 

H = -he^ 



G = -be"^ -4{H-G). (A.22) 



Note that H, H and G, G satisfy the same condition at t = — oo. 
On the other hand the derivatives d^H = h and d^G = g satisfy 



h = -he^g 

g = (4 - 6e^)5 - 4/1. (A.23) 

The conditions at t — > — oo are vanishing for both h, h and g, g. 
Introducing the energy £ = -^H^ + ^e we get: 

£ = 6ae^e2(*"i) < 

and hence 

6e^(*) < S{t) < S{-oo) = 2. (A.24) 

Therefore ^ > as far as /i < and /i < as far as ^f > 0. This conditions are indeed verified for 
t K, — c)o SO that they are true for all the time. Then G is increasing as well as 

Z = — I dte'^. 



Proof of Theorem [T^ 

To prove the concavity of F we will prove that, for any ai, 6i, 02, 62 : 

Let a = 2ii22j and b = ^~^ ^ . By the linearity of Fafi{uj) as a function of a,b at fixed oj, we get 
that 

F{a,b) = Fa^b{uJa,b) = -^Fa^^bi{uJa,b) + ^Fa^^bii^a^b) 

111 

> 2^a^,bj(^a^M) + 2^a2,b2i^a2,b2) = -(i^(ai, 6l) + -F(a2, 62)), 

where we used the fact that ujai,bi and U!a2,b2 are the minimizer for Fa^^^bi and Fa2,b2 respectively. 

The smoothness of F comes from the fact that the solution of the canonical variational principle 
depends smoothly upon a, b. By taking the derivative of F with respect to a we get 

dF dFa,b{(^a,b) J, X 

oa da 
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Here we have used the fact that the derivative of F with respect to u; evaluated in u)a,b vanishes, 
and the fact that the derivative of Fa^b with respect to the parameter a is given by /. In the same 
way we get dF{a,b)/db = —Eiuja^h)- 

Finally, the concavity of Fa^b implies that dl /da > 0, and that dE/dh > again with the 
notation I{a,b) = I{uja^b)-, E{a,b) = E{uJa,b)- 

Now it remains to prove iii). Again, the concavity of F{a, b) implies the existence of the convex 
function S*{I,E) defined in fO]) . Now we want to prove that S{I,E) = S*{I,E). First of all let 
us notice that 

S*{I, E) = sup (F{a, b) + bE + al) = F{d, b) + bE + al, 



where a, b is the unique maximum point for S*{I, E). Therefore, for any a, b 

S*il, E) > F{a, b) + bE + aI = S{uJa,b) - b{E{uJa,b) - E) - a{I{uJa,b) - I)- 

Now, since F is concave and smooth, we know that for any /, E, there exists unique a, b such that 
I(iJa,b) = Ii and E{uJa,b) = E. By choosing a, b in this way in the previous equation we get 

S*{I,E)>SiLOa,b)>S{I,E). (A.25) 

On the other hand, let ujj. : k = 1,2,... be a minimizing sequence for S{I,E), and u a limit 
point for it. By lower semicontinuity of S we know that S{uj) < S{I, E). Therefore, for any a, 6, 

S{I, E) > S{lo) = S{uj) - bE{uj) - al{uj) + bE{uj) + al{iu) 

> F{a, b) + bE{uj) + aI{oj) > F{a, b) + bE + al, (A.26) 

where we have used the continuity of E from which E{uj) = E, and the lower semicontinuity of /, 
from which I{io) < L 

Since a,b are arbitrary in ()A.26p . we get S{I,E) > S*{I,E),. Since we have already proven 
fOS]) this yields S{I,E) = S*{I,E). 

Finally let us notice that 

S{I,E)=S*{I,E)=S{u;-,-,), 

where a,b is the unique minimun point for S*{I,E), where the relation between a,b and I,E is 
smooth and bijective (equivalence of the ensembles). 
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